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Abstract 

We show in a rigorous way that Crum's result regarding the equal eigenvalue spectrum of 
Sturm-Liouville problems can be obtained iteratively by successive Darboux transformations. 
Furthermore, it can be shown that all neighbouring Darboux-transformed potentials of higher 
order, Uk and Uk+i, satisfy the condition of shape invariance provided the original potential u does 
so. Based on this result, we prove that under the condition of shape invariance, the n th iteration 
of the original Sturm-Liouville problem defined solely through the shape invariance is equal to the 
n th Q rum transformation. 
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INTRODUCTION 



Supersymmetric quantum mechanics 
transformation 



11 a, 



the factorization method 



the Darboux 
to- 



Ill 



4J, Crum's generalization of the former results [5|, the isospectral Hami 
nians based on the Gelfand-Levitan equation f], ?], Isl, Q] or the Marchenko equation 
and the shape invariance condition on the potentials [3] together with a transformation 
defined through this condition have been in the last two decade an active area of mathe- 



matical physics 



16 



17 ] and pure mathematics [18|,ll9|,|20|. The main concern of 



these areas has been the construction of isospectral Schrodinger operators and the analytical 
solvability of the Sturm-Liouville problem. The field allowed a deeper insight into the eigen- 



value problem and served as a source for many new ideas and generalizations 



21 



22, 



23 



24]. 



Indeed, it is almost impossible to quote all research papers on the subject (suffices to note 



that one review [25( and several books have been devoted to the subject 



26 



21, 



28 



29 



3Q|). 



The applications ran ge f rom constructing new solvable potentia 



differential equations 



3l|, atomic physics 32], nuclear physics 



s in quantum mechanics, 



33] , classical mechanics 34] , 

38|, 



acoustic spectral problems [35| to quantum gravitation [36j, |37J and neutrino oscillation 
to mention a few important areas. 

Mathematically, not all these transformations mentioned above are equal, or at least this 
is not apparent at first sight. For instance, the usual Darboux transformation is not the most 
general solution of the Riccati equation and as such does not give us the most general trans- 
formation in connection with the isospectral eigenvalue spectrum. On the other hand, the 
generalization of the Darboux transformations, namely, the so-called Crum transformation 
appears to be much more complicated than the original Darboux result and as such seems to 
offer us new avenues to construct new potentials. The third transformation of a Hamiltonian 
which we have in mind (defined here in equation (92)) is closely related to the condition of 
shape invariance. Hence, without doubt, there is some need to at least classify these trans- 
formations according to the complexity or generality and to uncover their relations between 
them. One such result in this direction is the nonequivalence of the Abraham-Moses [7] and 
Darboux constructions shown in j^]. Two remarks are in order here. Firstly, it is understood 
that unlike the Darboux transformation, any transformation in connection with the shape 



2 



invariance is, of course, limited to the set of shape invariant Hamiltonians. Secondly, for 
completeness it is worth noting that the level of complexity of isospectral quantum systems 



can 



39. 



3e increased by considering non-linear and higher order supersymmetric transformations 



40 



41 



42]. These are transformations which cannot be reached by iterative Darboux 



transformations. In this work, however, we will not consider these kind of transformations 
and restrict ourselves to the Darboux case. After some preparatory statements we will show 
that the undertaking to uncover relationships between the transformations gives a simple 
result, namely, allowing the use of higher order Darboux transformations, we can state that 
all three transforms of the original Sturm-Liouville problem are equal. This result is based 
on a theorem which we prove in the present paper concerning higher order Darboux transfor- 
mations of shape invariant potentials denoted by w D [/c]. The theorem states that provided 
the original potential satisfies the shape invariance conditions, all pairs u D [/c], u D [k + 1] are 
also mutually shape invariant. The theorem can be proved by induction. Interestingly, it 
intertwines this induction with another statement, this time for the wave functions. We 
illustrate the theorems by two examples. 



II. CRUM'S RESULT 

In this section, we briefly present Crum's result and comment on one identity on which 
Crum's result is partly based. This identity is crucial for the subsequent results which we 
will elaborate upon in the next section. 

Let 

W k = W(il> 1 ,ih,...,il> k ) = debA , ^ij = -^f i,j = l,2,-,k, (1) 
be the Wronskian determinant of the functions i/)i,ip2, ■■■,'4 ) k, an d 

W k , a = W(^i,ih,-,Ms) ■ (2) 

Theorem 1 [Crum]. Ifipx,^, ■■■,i>n are the solutions of the regular Sturm-Liouville prob- 
lem 

- + mp s = \s4> s (3) 
then ip c [n] satisfies the Sturm-Liouville equation 

- + u c [n] ^[n] s = X s ^[n] s , (4) 



with ip [n] s and u [n\ s given by 



ip s — ► ip C [n] s = 



W„ 



w n 



and 



u — > u c [n] 



d 2 

u-2—\nW n 

dx z 



(5) 



(6) 



Note that the Crum transforms of ip an d u are not defined iteratively. By [C] we wish 
to distinguish the Crum transformation from other transforms (like Darboux) which will be 



and 28]. We 
5| which we will 



defined later in the text. The proof of Crum's theorem can be found in 
comment here only on one cornerstone of the original proof given by Crum 
also use later. The first step in the proof of Crum's result on the Wronskian determinant is 
to consider the derivative of W/~. Taking the derivative of 



Wo 



dx dx 



(7) 



we find the rather obvious result, 

dW 2 

dx 

where we used the notation tp" for ^Jr-, 
the n x n case. 



(8) 



^1 ^2 

1,2. This result can be readily generalized for 



Lemma II. 1. For the derivative of a Wronskian determinant we have 

W > n = {^Mji, + 4 n) Mg } + ■ ■ ■ + M) + #<n) } ■ 0) 

Assume the result to be valid for n — 1. Using the Laplace expansion according to the 
last line of the Wronskian W n we get 



W' 



#- X) «o)' + " " " + ^t? « 1>n) )' + ^ «„))' 



(10) 



where every (n — 1) x (n — 1) determinant M^K is a Wronskian for which, by assumption, 
the theorem is valid. Hence 



is a determinant whose two last lines are equal and therefore detB = 0. The result (Tl0|) can 
be written as 



w' 



ft ft 2 



ftr 2) ^r 2) 



4>n-l 4>n 



in) 



4«» 



(12) 



We can now state a result which will be of some importance later and which is one of the 
important ingredients in proving Theorem 1 of Crum. 

Lemma 1 [Crum]. The Wronski determinant of the two Wronskians, W n and W n -i jS , is 
equal to W n , s W n -x. In other words 

W(W n} W n . 1 , s ) = W ns W n . 1 . (13) 

The proof relies on the Jacobi theorem for determinants (see Appendix A). We refer the 
reader to Appendix A for the proof of this Lemma too. 

It is well known that for n — 1 the Crum transformations reduce to the Darboux trans- 
formation when Wi = ft, and W\ tS = W (ft, ft). Specifically, we have 



<ni] s =v c [i] 



Wl )S , ft 

—— =V>s- — 8 > 1 , 
ft ft 

u D [l] = u c [1] = u - In W 1 = u - ■ 

ax 1 ax ft 

We can define higher order Darboux transformations iteratively by 
u d [k _ 1] ^ u D [k]=u[k _ 1] _ 2 d l k - !]ft) 



(14) 
(15) 



tP D [k - 1} S i, D [k] 



D 



dx tp D [k-l] k ' 
Wty D \k-l] k ,^ D [k-l] s ) 



s > k 



(16) 



r [k - i] k 

Obviously, the last equation can be written also in a way which resembles more the first 
Darboux transformation, i.e., 



r [k] s = ^ [k - i] s y 



(r 5 [k - i] h y. 



r [k-i\ t 



(17) 



r [k - i} k 

It is a priori not clear as to what connection the k-th Darboux transformation has with 
the fc-th Crum transformation and if they can be related at all, except for the definition at 
the lowest order of Crum's transformation. The answer is provided in the next section. 



III. THE CONNECTION BETWEEN HIGHER ORDER DARBOUX AND CRUM 
TRANSFORMATION 



To this end, let us first examine the simplest case of k — 2: 



u D [2] = u D [1] - 2 



dx i> D \l\ 



(18) 



Since w[-D][i] is the Darboux transformed potential the above equation reads 



u D \2\=u-2^-[?±+ vr n L r T 7 ) . (19) 

[ J dx Ui ^ D [i] 2 ; 1 ; 

According to (19), V'" [1] 2 — an d on account of the simple identity W n ^ n+ i = W n+1 , we 
can write, 



u D [2] = u - 2 



which finally gives 
Similarly, the eigenfunctions 



(20) 



(21) 



^ [2] s = 



^ [1] 2 ^ [l] s 

ib D [i] 2 [i]. 



^ [i] s 



(22) 



can be cast into the form 





w 2 


w 1>8 


J_ 

















w 2 w. 



W 2 \ I Wi, s 



w 2 



Wo 



(23) 



With the help of the standard property of determinants, namely det (z±, z iy z n ) 
det (zi, Zi + azk, z n ) the last equation reduces to 



r [2], 





w 2 


w 1;S 


J_ 












W{ 





Wo 



(24) 
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Applying the result of Lemma 1 and remembering that W\ = ipi , one finally finds 



r [2], 



Wo 



(25) 



The steps above will serve as a beginning of the induction proof of the following general 
statement: 

Theorem III.l. The n-th Crura transformation is equivalent to the n-th higher order Dar- 
boux transformation. This is to say, any Crum transformation can be reached iteratively by 
successive Darboux transformations, i.e., 



u c [n] = u D [n] , 
i> C [n] s = ip D M s • 

Proof. Assuming the theorem to be valid for n means that the statement 

d ^ D [n - 1]J 



u D [n] = u D [n - 1] - 2 



dx ip D [n — l\ n 



^ [n] t 



W{^ D [n-l] n jn-l] s ) 
V> [n _ i] n 



s > n 



is equivalent to 



/ D r l W n s 
W [n\ s 



s > n 



Based on that, we have to show 

u c [n + 1] = u D [n + 1] = u D [n] - 2 

and 



d 0> D [n] n+1 Y _ d W' n+1 

— xt — 



dx ip D [n] 



n+l 



dxW n+ i 



The validity of the hypothesis of the induction for n allows us to write 

d [n] n+1 ) ! 



u [n + 1] = u [n] - 2 



dx ip D [n 



n+l 



dx \ W, 



+i 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 
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Using the validity of the hypothesis for n, but this time for the wave functions, implies 

/ . fw n+ i\'\ 



u D [n + 1] = u - 2 



d_ 

dx 



w; \ w < 

— Vl + _i 

Wn 



\ 



u 



I 



Similarly, the result for the eigenf unctions may be written as 



w, 



n + l 



It', 



n+l 



w n 



ip D [n] 



w, 



n+l 



n + l 



One easily proceeds now to verify the validity of the following equation 



i> D [n + l] s = 



By virtue of Lemma 1 we can assure that 



w n+1 w n , s 


1 

w n 


w n+l W n , s 


W^+i w^w n+1 w„' iS w n w n ,s 

W n Wl W n Wl 


w; +1 w>, s 


w n+1 




w n+1 



(33) 



(34) 



(35) 



n+l 

is true which completes the proof. 

It is instructive to follow this theorem by an explicit examples. 

IV. TWO EXAMPLES 

In this section we will demonstrate the above theorems by two examples. We choose first 
a potential which satisfies the condition of shape invariance (Morse potential) followed by a 
simple example which falls into the class of non-shape invariant, but solvable potentials. 

Let us consider, as an example, the Sturm Liouville problem with the Morse potential, 
i.e., 



A 2 - A ( A + ) sech 2 (n.r) 



u (x; A) = 2 

The super-partner of this potential corresponding to ipi and Ai is 

u c [1] = u D [1] = 2 [A 2 - A4 lS ech 2 (ax)] 
and the first three eigenfunctions are given by 



(36) 



(37) 



V2A 

1. ipi = ci [sech (ax)] a 



2. ip2 = C2 sinh (ax) ipi 



3. ^3 — C3 ^— cosh 2 (ax) + ( 2v/2 "^ a ) ginh 2 («x)j ipi 



In the following we will not determine the constants Cj as they are of minor importance 
for our results. Secondly, the results become increasingly complicated. For instance, to 
calculate c\ we can use 

i , , 1 „ I V2A 1 1 a \ , 00 . 

sech(ax) « ax = =— 2 -ri , -, H , r- (3o) 

1 ; v/2A \ a '2' a ' v^4/ 1 ; 

where 2-^1 is the hypergeometric function. The corresponding eigenvalues can be compactly 
written as 

K = 2{a--(a-^-^\\ . (39) 



It is convenient to define A„ as 



y/2 



not , . 

A n = A — — (40) 



such that the eigenvalues read now 

X n = 2(A 2 -A 2 n _ 1 ) n = 1,2,3,... . (41) 

The first three are explicitly given as follows 

Ai = , A 2 = 2v / 2Aa - a 2 , A 3 = 4V2Aa - 4a 2 . (42) 

Besides equation (|37j) we will also need the following results: 

t/> c [1] 2 = ip D [1] 2 = c 2 a cosh (ax) V>i, (43) 
^ c [i] 3 = t/, d [l] 3 = 4^03^! sinh (ax) cosh (ax) ^1 . (44) 

To show explicitly the equality u c [2] = u D [2], we start with u c [2], i.e., 

M G [2] = u-2— lnW 2 . (45) 
ax 2 

Since Wj 2 = W 2 and using 

' = c 2 a cosh (ax) ^1 (46) 



we have, 



but also 



In W2 = In C2OL + In cosh (ax) + 2 In 
d 2 



— In W 2 = a (a - 2\[2A^ sech 2 (ax) . 
Finally, with fl36|) we arrive at 



u 



c 



[2] = 2 [A 2 - AxAasech 2 ( 



ax 



Next we turn to the expression for u D [2], namely, 
Taking into account equation ( l4~3l) we obtain 



(JjD hi y 

vy n L r j2; = -y/2Ai tanh (ax) 



and therefore 



dx i) D [l] 2 

From this we conclude (see (1501) ) that 



•v^aAisech 2 



[2] = 2[A 2 - Assert! 2 



[ax 



(47) 
(48) 

(49) 

(50) 

(51) 
(52) 

(53) 



and hence 

u c [2] = u D [2] . (54) 

The transformed potentials here have almost identical functional form. This is, of course, 
due to the choice of the potential and need not be so in other cases. 
To demonstrate that ip c [2} s = ip D [2} s , we calculate ip c [2] 3 to be 



^ [2]. 



Wo 





"02 


^3 


-Ai^i 




03 










^ i>2 





^2-rrr- ~ A 3 3 • 

v<2 



(55) 



Making use of 



A 2 ^: 



W h3 /V2A2C3AA . 



W 2 



a 



sinh (ai) ^1 



(56) 
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this becomes 



On the other hand 



tp c [2] 3 = A3C3 cosh 2 (ax) ip\ (x) 



(57) 



4> D \% 



^[l] 2 ^[i] 3 

(vw (m) 



(^[1] 2 ) 2 



A3C3 cosh 2 (ax) ^1 



v-[i] s 



V»[i] s 



(58) 



where on the right hand side we already dropped the distinction between D and C (see (JI2 
and (j44p ). The simple conclusion that we can draw is 



i, c [2] 3 = ^ [2] 



(59) 



It is instructive to consider also a case of a solvable, but non-shape invariant potential. 



Many such cases are known (see 



43 



44 



45] and the discussion in [25j) and explicit proofs 



that these potentials fail to satisfy the shape invariance condition were given. For instance, 



for the case of the Natanzon potential this was shown in [4a]. Many of these potentials are 
complicated and some, like the Natanzon case, only known in implicit form. Therefore, for 
the sake of efficient calculations, it is recommendable to develop first a fast algorithm to 
perform the desired calculations. We will do exactly that before explicitly giving the explicit 
example of the Ginocchio case. Imagine we would like to calculate ip D [2] 3 . In turns out that 
the calculation can be greatly simplified by invoking the ratios h n = ip' n /ip n where ip n is, as 
usual, the eigenfunction to the e n eigenvalue. It is now a straightforward exercise to show 
that 

(^[iy (nm 



r [2] 



is equal to 



(hi - h )ipi 



(hi - ho) {(h 2 - h )' + (h 2 - ho) h 2 } ^1^2 
(h 2 - h ) {(hi - h Q )' + (h - h ) h} ^1^2 



Using the Schrodinger equation the latter simplifies to 

[e {h 2 - hi] - ei {h 2 - h } + e 2 {hi - h }] 



(hi - ho) 



i>2 



(60) 

(61) 
(62) 

(63) 
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It is obvious that, provided we know the functions ho, hi, hi and ip2, this expression allows a 
fast calculation or the wave function ■0 £> [2]3 for arbitrary potential. The Crum's result gives 



w. 



2..°, 



Wo 



Wo 



(/ii - ho)ipotpi 



(64) 



where W-* is 



W, 



v>0 


4>1 


^2 








^2 




^0 














h 4>o 


hiipi 


h 2 ip2 






hii>i 


h 2 ip2 






1% 




< 








u - e ^o 


u - e^i 


u - ei-02 



4>0 V>1 ^2 
/i ^0 ^lV'l ^2^2 

e ^ ei^i £2^2 



(65) 



Hence, taking ( 1641) into account, we can show that 

,cr 9l [go {^2 - M - ei {fr 2 - M + e 2 ~ M] ^1^2 , , 

^ [Js " (/n-fc>)tMb ' ^ j 

which obviously implies that = ^ [2] 3 . This, as it stands, is a general proof for a 

sub-case of our general theorem. On purpose above we have used different steps than in 
the proof of our general theorem. The idea behind it is to demonstrate that in an explicit 
example we would be only repeating he very same steps as above. It is therefore sufficient 
to calculate every time only the right hand side of ( 1661) . The equality ^ D [2]3 = ip [2)3 
is guaranteed by ( |63j) and [66]). We can now apply the results for V' D [2]3 by choosing the 
Ginocchio potential 



V{x) 



{ -fv (v + 1) + I (1 - /3 2 ) [5 (1 - /3 2 ) y* - (7 - /? 2 ) y 2 + 2] } (l - y 2 ) 



(67) 



where y (x) satisfies the following differential equation 



(68) 



and (3, v are parameters. 

The wave functions of this problem are known to be expressible through Gegenbauer 
polynomials Cn\x), namely 



^ = (1 - [9 (y)]- (w)/4 (/ (y)) 



(69) 
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where 

g(y) = i-(i-P 2 )y 2 , f{y) = ^M= (70) 

V9 (y) 

The value of /i n is connected to the eigenvalue e n by e n = — /U 2 /3 4 such that 

/in/3 2 = v / /5 2 (^ + l/2) 2 + (l-/3 2 )(n + l/2) 2 - (n + 1/2) . (71) 
The first four Gegenbauer polynomials are given as follows 





f-1/2) 


(/(v)) 




H/2) 


(/(y)) 




f-1/2) 


Cf(y)) 


£,0*3+1/2) 


(/(v)) 



(72) 



| (// 3 + 1/2) (/i 3 + 3/2) (Aia + 5/2) [/ (y)] 3 
-2 (/i 3 + 1/2) (/i 3 + 3/2) [/ (y)] . (73) 



These functions can be used, in the next step, to compute explicitly the ratios hi = 4>'i/ipi. 
We obtain 

ipo g{y) 

hi = ^ = kM + ^ = (I^!> {1 _ 2(1 _ /3V} 

V»i \g{y)\ [f{y)\ y 

h2 = ^ = Wl 2 [/(,)][/(,)]' 

* V s ([/fo)] 2 -<^j 

where we have used 

VM = v {1 - v) - (75) 

Noting that the hi are proportional (1 — y 2 ) and that h\ — h = (1 —y 2 )/y we can insert our 
results into equation fl66l) to obtain 

V>°[2] 8 = (l-/3 2 )^ /2 ( / , 2 + l/2)(2 / , 2 + 3)[^(y)]-^ +1 )/ 4 

x { (c 2 - e ) ([f (y)] 2 - {2 ^ +3) ) - (ci - ^o)2[/(i/)] 2 } • (76) 

Turning our attention to the potential the superpartner of V in equation fl67|) it is not 
difficult to see that the superpartner is given by 

V D [1] = V c [l] = V-2^lnW 1 = V-2P-^f 

(lx- dr g (//) 

= V + 4 (1 - (3 2 ) (1 - 3y 2 ) 2 (1 - y 2 ) [1 - (1 - (3 2 ) y 2 ] . (77) 
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The second Crum iteration yields 



d 2 



= V 



{-2+ (I-/? 2 ) [5y 2 - 3] } + 10y 2 (1 - (3 2 ) 



(l-y 2 ) [l-(l-(3 2 )y 2 ] (78) 



To proof that this is equivalent to the second Darboux transformations it is convenient, as 
it was the case with the wave functions, to provide first a general proof for this sub-case. 
Starting with the definition, it is straightforward to show that 



^ D [1] 2 = (/H- Wi, 



(79) 



which leads to 



V° [2] = V D [1] - 2^^P = V - 2^ In ^ - 2-£ In ( hl - h ) 



dr 



dr 2 



{In (h - h ) ^i} = V - 2^ In W 2 = V c [2] 



(80) 



In taking explicit examples we would be only repeating the very same steps as above. This 
demonstration concludes our two examples. 



V. THE CONNECTION BETWEEN SHAPE INVARIANCE AND CRUM 
TRANSFORMATIONS 

In view of the results of the previous section we can now drop the distinction between 
higher order Darboux (D) and Crum (C) transformations. 

Let a denote a set of parameters in the original potential, i.e., 

u = u(x; a) . (81) 

The condition for shape invariance of u is given by 

u[l}(x;a)=u(x;f(a)) + R(f(a)) , (82) 

where u [1] {x\ a) is the first Darboux transform of the original potential, / transforms a into 
another set / (a) and R (/ (a)) is a function of the parameters. In the following, we use the 
usual notation a m = f m (a), where m indicates the function / applied m times. 
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In the preceding section we established an equivalence between higher order Darboux 
transformation and the Crum result. Since the shape invariance is given in terms of the first 
order Darboux transformation, it is legitimate to ask if higher order Darboux transforma- 
tions (Crum transformations) play a role in the Schrodinger equation with shape invariant 
potentials. As a first step we will prove the following theorem 

Lemma V.l. Under the condition of shape invariance one has 

ij) a (ar; oi) = if) [l] s+1 (ar; a) . (83) 
up to a multiplicative constant and 

X s (a 1 ) + R(a 1 ) = X s+1 (a) . (84) 

In the above ip s (x\ a) denotes the eigenfunction to the Hamiltonian with the potential u(x; a) 
with the eigenvalue X s . 

These results are not new. But since we will make use of them, we offer here a short 
proof. We start with initial Sturm-Liouville problem 

d 2 



dx 2 



+ u (x; a) I ip s (x; a) = X s (a) ijj s (x; a) (85) 



and 

d 2 



+ u[l](x;a))^[l] s (x;a) = X s (a)^[l] s (x;a) s>l. (86) 



dx 2 

Equation ( 1851) is valid for any a, hence we may write 

d 2 



+ u (x; ax) %l) s (x; ai) = X s (at) ip s (x; a{) (87) 



dx 2 

and add R (f (a)) ip s (x\ f (a)) on both sides implying the following identity 
d 2 



dx 



2 



+ u (x; ai) + R (ai) ip s (x; a%) = (X s (ai) + R (ai)) ip s ( x ', a i) ■ (88) 



Due to the shape invariance condition (for the sake of formulating the next Lemma we can 
say that u [1] and u are pairwise shape invariant) this becomes 

d 2 



^ 2 + u [1] (x; a) ) ip s (x; at) = (X s (at) + R (at))ip s (x; a x ) ■ (89) 
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Without loss of generality, the spectrum can be ordered as Ai < A2 < A3 < .... Hence, 
{A s (a)}, {A s (ai)} and {A s (ai) + R(a\)} are similarly ordered sets, ip [l] s+ i is then an eigen- 
function to the ordered spectrum A 2 < A 3 < .... We can conclude that up to a multiplicative 
factor 

ip s (x;ai) = i) [l] s+1 {x, a) (90) 

and 

X s (a 1 ) + R(a 1 ) = X s+1 (a) . (91) 
In preparation of the main theorem of this section we prove the next Lemma. 

Lemma V.2. By virtue of the the above Lemma and under the condition that u and u [1] 
are pairwise shape invariant, u[l], u [2] are also pairwise shape invariant i.e. 

u[2](x;a)=u[l](x;a 1 ) + R(a 1 ) . (92) 

The proof can be done in two steps. 

1. The condition of shape invariance and the definition of the Darboux transformation 
allows us to write 

u (x; a) - 2 ( g >°) = u hi f X] a \ = U U ; a A + r f a A (93) 
ax ipi (x; a) 

which remains valid if we replace a by ax, i.e., 

u (x; ax) - 2- — — = u 1 fx; a{\ = u(x; a 2 ) + R (a 2 ) . (94) 

ax ipx {x; ax) 

Hence, we easily obtain 

/ \ r n / \ d ib'i (x: ax ) /„_ s 
„( I;ai )=«[l](x; 0l )+2 s ^^. (95) 

2. Applying the Darboux transformation ffl6|) . once again on u [1] (x; a), gives 

W 2] M = M1] M -2^M. (96) 

On the other hand, using the shape invariance condition leads to 

u [2] (x; a) = « (1; a x ) + R (ax) - 2 ± ^^^ a) . (97) 

ax ip [1\ 2 [x; a) 
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The result in the first step of the proof, (195]) . can be used to derive the following 
equation: 

u [2] (x; a) = u [1] (x; ax) + ( fll ) + 2 A ( - ^ ( * ; Q . } 1 . (98) 

ctx ^i(x;ai) ■?/> [1J 2 (x; a) J 

If we now apply (190]) from Lemma V.l for s — 1, i.e., 

^ [1] 2 (x; a) = ^! (x; a x ) (99) 
we obtain the desired final expression which we wanted to prove, namely, 

u [2] (x; a) = m [1] (x; ai) + (ai) . (100) 

For the sake of a more compact notation of the properties of the potential and wave 
functions, let us now call the property fl90|) shape invariance for eigenfunctions (or better 
the two eigenfunctions involved are pairwise shape invariant) and shape invariance for 
the eigenvalues. Note that the shape invariance of the wave functions follows from the shape 
invariance of the potentials. From the shape invariance of the eigenfunction we can, in turn, 
conclude that the next two pairs of Darboux transformations of the potential are pairwise 
shape invariant. One is led to the conjecture that the chain continues: from Lemma V.2 
one can show that the next pair of higher order Darboux transformations of eigenfunctions 
are also pairwise shape invariant, from which it follows that the next higher order pair of 
Darboux transformed potentials is also pairwise shape invariant. Indeed, we can prove the 
following theorem extending hereby the notion of shape invariance. 

Theorem V.3. All neigbouring higher order Darboux transformed potentials and eigenfunc- 
tions are pairwise shape invariant. This is to say, 

u [k] (x; a) = u [k - 1] (x; a x ) + R [ax)) (101) 

and 

ij[k} s+1 (x;a)=ij[k-l} s (x;ax) , (102) 

up to a multiplicative factor. In more detail, $101\) implies which, in turn, implies 

u [k + 1] (x; a) = u [k] (x; ax) + R (ax) ■ (103) 
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The proof proceeds via induction whose first step consists in Lemma V.2 and Lemma V.l 
or in fl90l 11001) . We assume the hypothesis of the induction to be [ fjlOXp =>- (11021) ]. This is 
sufficient since we start with the original shape invariance condition for potentials and the 
first step of induction is presented in Lemma V.l and Lemma V.2. We have to show that 
under this condition 

1>[k + l] s+ i 0; a) = -0 [k] a (x; ax) (104) 

holds, from which, in turn, 



u [k + 2] (x; a) = u [k + 1] (x; ai) + R (ai) (105) 



follows. 

1. We have 

d 2 \ 

— + u [k] {x; a) j ip [k] s {x; a) = X s (a) if) [k] s (x; a) s>k. (106) 

2. Since the above equation is valid for any a, it is also valid when a is replaced by f(a). 
If we add R (ai) ip [k] s (x; a{) on both sides and make use of the induction hypothesis 
we arrive, for s > k, at 

d 2 



(J )2 +u[k+l] (x; a) J V [k) s (x; ai) = (A s (ai) + R (ai)) ^ [A;], (x; a x ) . (107) 
Equation then tells us that, 

1>[k] s {x; ai ) = ip[k + l] t+1 {x;a) , (108) 
up to a multiplicative factor. 
3. By definition we have 

u[k + l] (x; a)=u [k] (x; a) - 2^- ( -. ^+^ X]a A (i 9) 



dx\ip [k] k+1 (x- a) 
for any a. Hence also: 

„[» + i] („ ai) = » w (x; 0l) - 2 1 ) . (no) 

Taking w [&] (x; ai) from this equation and inserting the result in the induction hy- 
pothesis, one can easily show that 

d fip'[k] fc+i [ x 'i a i) 



»[* + !] (*; a) = u [k + 1] („ a,) + 2^ ^ — fe ' | + R (a.) . (Ill) 
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k+2 y x 'i a ) 



4. Again per definition we know that 

u [k + 2] (:r; a) = u [k + 1] (x; a) - 2-f- I r ; r L ,~ ■ f^;^ ) . (H2) 
L J V ; 1 J V 1 dx V ^[& + l] fc+2 (a;;a) ' 

5. Combining the last two equations yields 
u [k + 2J (x; a) = w [« + 1J (x; ai)+R [a 1 )+2— 1 



dx \ip [k] k+1 (x; ax) ip[k + l] k+2 (x; a) 

(H3) 

6. The last step consists in using the already established result (I108[) to obtain 

u [k + 2] (x; a) = u [k + 1] (ac; ai) + i? (a x ) , (114) 

which completes the proof. 

The shape invariance condition (more accurately, the shape invariance between u and 
u[l]) allows one to define a new Hamiltonian of the order s: 

d 2 s 
H^ = -—+u(x;a s )+J2 R ( a k) ■ (115) 

k=l 

Note that this definition makes no reference to higher order Darboux (or Crum) transfor- 
mations. However, by virtue of the Theorem V.3 we can iterate 

d 2 5-1 
Hi 1 = -— + u[l}(x;a s ^) + J2 R ( a k) 

k=l 

d 2 s ~ 2 
= ~~^2 + «[2](x;a s _ 2 ) + ^2R(a k ) 

k=l 

= ... = --^ + u[s}(x;a) = H? . (116) 
Hence, in view of the above and the Theorem III. 1 we can state as a corollary 

Collary V.4. Under the condition of shape invariance all three transformations are equal, 
i.e., 

H SI = H D = H C_ (n7) 
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VI. EXAMPLE 



We will continue here with our example of the Morse potential, now however emphasizing 
the aspect of shape invariance around Lemmas V.1-V.2 and Theorem V.3. Indeed, the Morse 
potential is shape invariant. One defines the action of / as 

f(A)=A x = A-^= (118) 

in accordance with the notation in (HOT) . R is identified with 

R (A 1 ) = 2 (A 2 - A*) . (119) 

Note that 

fa (x; A\) = c\ (Ai) [sech (ax)] a = c\ (At) [sech (ax)] a 

= c cosh (ax) fa (x;A) (120) 

immediately leads to 

iP[l] 2 (x;A)=fa(x;A 1 ) (121) 
which is valid up to a multiplicative constant. One also verifies that the equality below 

fa (x; A x ) = c 2 (Ax) sinh (ax) fa (x; A x ) (122) 

= c 2 (Ai) sinh (ax) cosh (ax) fa (x; A) (123) 

together with (jUJ) has as a consequence the following identity (again up to a constant 
multiplicative value) 

i>[l} 3 (x;A) = fa(x;A 1 ) . (124) 

Equations (I12ip and (11241) are explicit examples of the result (183]) in Lemma V.l. Regarding 
the eigenvalues, i.e. the property ( |84l) in the same Lemma, let us first note that another 
compact notation for equation (14ip is 

A n (A,) =2(A\- Al) (125) 

which leaves us with the identity 

\ n (A 1 ) + R(A 1 ) = \ n+1 (A) , (126) 
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as it should be according to Lemma V.l. Finally, we can also give explicit examples regarding 
Theorem V.l. Due to the results from section 5, we can write, 

a \ 



u[l] (x; A 1 ) + R(A 1 ) 



Af- A 1 [ Ai- — J seclr (ax) 



+ 2 (A 2 - Af) 



2 [A 2 - A^sech 2 (ax)] = u [2] (x; A) . 



(127) 



This demonstrates in an explicit example the result fllOip from Theorem V.l. Last but not 
least, one sees that equation (1431) can be written as 



if) [1] 2 (x; Ai) = c 2 (A^ a cosh (ax) if)\ (x; A\) 



(128) 



which, according to (1121 j) can be cast into the following form: 

if) [1] 2 (x; Ax) = - 2 — — —a cosh (ax) if) [1] 2 (x; A) = — ^ 2 - — -a 2 cosh 2 (ax) if)\ (x; A) 

= ciP[2} 3 (x;A) (129) 

with c a constant. To arrive at the last result we have used (1571) from which one can also 
determine the constant c in terms of A 3 , c 2 (A), c 2 (Ai) and c 3 (A). Obviously, the above 
equation falls into the category of explicit examples of (11021) . Note that in none of the above 
examples we have used the actual Lemmas or Theorems to be exemplified (as it should be 
if an example carries some meaning). 



VII. CONCLUSIONS 

In the present work, we have clarified the relations between the Darboux and Crum 
transformations. We have shown that the latter can be reached iteratively by higher order 
Darboux transformations. This is valid for the potential as well as the eigenf unctions. If we 
subject the potential to the condition of shape invariance, another transform (not making 
use of Crum transformation for n > 1) is possible (equation (I115P ). We prove that this is 
also equivalent to the Crum transform. The main steps of this proof involved establishing 
( 11011) . (I102p and (18"1|) . The first result, namely (11011) . is a generalization of the original shape 
invariance condition. Note that ( 11021) and ( |84l) could be called shape invariance for the wave 
functions and eigenvalues. The results of the paper help to understand the relation between 
the different transformations of the Hamiltonian operator. Indeed, in view of our results, 
one could say that the Crum transformation which appears much more complex than the 
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original Darboux result is essentially an iterative higher order Darboux transformation and 
therefore not more complex than the former. 



VIII. APPENDIX A: AN APPLICATION OF JACOBI THEOREM 



The identity W (W n , W n -i jS ) = W ns W n -\ has been used by Crum in the proof of his 
theorem. We have also made use of it several times in the present paper. It therefore makes 
sense to sketch a proof of the same. 

Let us first establish some notations and definitions. Let A = [a i3 ] beanxn matrix. The 
determinant of A will be denoted by \A\ as usual. We call the minor M r , the determinant 
obtained by retaining from A the r lines ii, i 2 , i r and the r columns k±, k 2 , k r . One 
defines the complement of the minor M r as the determinant obtained from A by dropping 
the r lines ii, i 2 , i r and the r columns k±, k 2 , k r . This complement will be denoted 
by M r c . One then defines 



M< r > = (-1) 



ii+i2+---+i r +ki+k 2 +...+k r 



(130) 



Furthermore, let A be the matrix of the cofactors of A : 



A u A 12 
A 2 \ A 22 



A 2n 

A 



131) 



A n \ A n2 

and M r , M' r the minors of \A\ and A, respectively. 

Theorem [Jacobi]. With these notations, the theorem of Jacobi asserts that 

ML = |Ar^M (r) . 



(132) 



Before proceeding we make a small diversion to an example of the application of the 
above theorem starting with a Wronskian composed of ipi, tp 2 , ^3, ip s , i.e., 



\A\ = W : 



3,s 



Ipl i>2 ^3 

f\ ^2 ^3 IP's 

w m r s € 

< ^ ^ r: 



(133) 
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The matrix of the cofactors is then given by 



ti 2 tis tis 

til til til 
til' til' til' 



tii tis ti s 

til til til 
tit ^" €' 



+ 



^2 

< ^2 € 
K W C 



^2 ^3 
til V>2 ^3 
K < ^3" 



Y 2 


ibi 
Yi 


ib, 

Ys 




Yl 


ibi 
YS 


ib, 

Ys 




rl 


ibo 

Y2 


Ys 




Yl 


ibo 

Y 2 


tbi 

YS 


til 




til 


+ 


til 


til 


til 




til 


til 


til 


+ 


til 


til 


tis' 


til' 




til' 




til' 


til' 


til' 




til' 


til' 


til' 




til' 


til' 


til' 


ti 


ti 


ti 




ti 


ti 


ti 


+ 


ti 


ti 


i>s 




ti 


ti 


ti 


ti 2 


^3 
^3" 


tis 

til' 




til 
til' 


tis 

til' 


tis 
til' 




til 
til' 


ti 2 
til' 


tis 

til' 




til 
til' 


ti 2 
til' 


tis 

til' 


^2 


^3 


*ps 




ti 


ti 


^s 




ti 


ti 






ti 


ti 


ti 


^2 


^3 


tis 


+ 


til 


tis 


tis 




til 


ti 2 


tis 


+ 


til 


ti 2 


tis 


til 


til 


til 




til 


til 


til 




til 


til 


til 




til 


til 


til 



We choose as lines and columns: (11,12) = (3,4) = (ki,k 2 ) . Applying the Jacobi theorem 
gives us 



ti ti i> s 

til ti 2 tis 

til til' til 



ti ti ti 

til ti 2 tis 
til' til' til' 



ti 


ti 


*Ps 


+ 


ti 


ti 


ti 


til 
til 


ti 2 
til 


tis 

til 




til 
til 


ti 2 
til 


tis 

til 



W- 



3,s 



ti ti 
til ti 2 



(134) 



Using Lemma II. 1 the left hand side takes the form 



Aw? c — Wt 

dx z >* dx - 3 



W- 



2,s 



(135) 
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such that we can write 



dx yV2 ' s dx VV3 



-W 2 , a W 3 
Explicitly, this implies the following equality 



W 3 , a W 2 . 



(136) 



W 3 , S W 2 = W 3 -^W 2 ,, - W 2 , S -^W 3 = W (W 3 , W 2 ,,) 



(137) 



The proof of the general case does not require any new procedure and follows essentially 
the steps outlined in the example. Let W ns be the Wronskian of the n + 1 functions 
V>i, ■ ■ • ,ip n , V^namely, 



\A\ = W n , s = 



(138) 



^ 4 n) ... ^ 4 n) 

and A the matrix of the cofactors of W n:S . We would like to apply the Jacobi theorem for 
the choice 

(t u t 2 ) = (n,n+l) = {k u k 2 ) (139) 
such that r = 2. In this case we need 



ML 



= W(W n ,W n . hs ) , 



(140) 



-w n _ 1)a +w n 

where we have used explicitly the result of Lemma II. 1. Clearly, we have, 

M« = W n . x , 

such that the Jacobi theorem for the Wronskian A can be stated as 

W(W n ,W n - 1 ,,)=W na W n - 1 , 



(141) 



(142) 



which proves Lemma 1. 
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